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ABSTRACT 

We have undertaken a deep (c ~ 1.1 mJy) 1.1-mm survey of the z = 0.54 cluster 
MS 0451.6—0305 using the AzTEC camera on the James Clerk Maxwell Telescope. We de- 
\ tect 36 sources with S/N^ 3.5 in the central 0. 10 deg 2 and present the AzTEC map, catalogue 

and number counts. We identify counterparts to 18 sources (50%) using radio, mid-infrared, 
Spitzer IRAC and Submillimeter Array data. Optical, near- and mid-infrared spectral energy 
distributions are compiled for the 14 of these galaxies with detectable counterparts, which are 
expected to contain all likely cluster members. We then use photometric redshifts and colour 
| selection to separate background galaxies from potential cluster members and test the reli- 

ability of this technique using archival observations of submillimetre galaxies. We find two 
potential MS 045 1 —03 members, which, if they are both cluster galaxies have a total Star- 
ke \ formation rate (SFR) of ~ 100 MQyr -1 - a significant fraction of the combined SFR of all 
■ the other galaxies in MS 045 1 —03. We also examine the stacked rest-frame mid-infrared, mil- 
limetre and radio emission of cluster members below our AzTEC detection limit and find that 
the SFRs of mid-IR selected galaxies in the cluster and redshift-matched field populations are 
comparable. In contrast, the average SFR of the morphologically classified late-type cluster 
population is ~ 3 times less than the corresponding redshift-matched field galaxies. This sug- 
gests that these galaxies may be in the process of being transformed on the red-sequence by the 
cluster environment. Our survey demonstrates that although the environment of MS 045 1 —03 
appears to suppress star-formation in late-type galaxies, it can support active, dust-obscured 
mid-IR galaxies and potentially millimetre-detected LIRGs. 

Key words: submillimetre - galaxies: clusters: individual (MS 045 1 —03). 



1 INTRODUCTION 

Galaxy clusters are highly biased environments in which galax- 
ies potentially evolve more rapidly than in the field. The galaxy 
* E-mail: j.l.wardlow@durham.ac.uk populations of local massive clusters contain mainly early- 
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type galaxies which define a colour-magnitude relation (CMR) 
dVisvanathan & Sandage|[l977l : lBower et al J 1992b . However, stud- 
ies of clusters out to z ~ 1 suggest that they contain increasing 
activity at higher redshifts due to a growin g fraction of blue, star- 
forming galaxies dButcher & OemleJ 19841) . Over the same redshift 
range there appears to be a growing deficit in the CMR popula- 
tion at faint magnitudes, as well as a claimed increasing decline in 
the numbers of SO galaxies, suggesting that the blue, star-forming 
galaxi es may be transforming into these passive populations with 
time dPressler et alj|l997l: ISmail et alj|l998l : be Lucia et ai]|2007l ; 
IStott et al.ll2007l : lHolden et al.ll2009h . 

The blue, star-forming populations within the clusters are ac- 
creted from the surrounding field as the clusters grow via gravi- 
tational collapse. The evolution in the star-forming fraction in the 
clusters may thus simply track the increasing activity in the field 
population at higher redshifts. The increasing activity in the field is 
also reflected in an increasing number of t he most luminous (an d 
dusty) star-burst galaxies with redshift (e.g. lLe Floc'h et al]|2005l) : 
the Luminous Infrared Galaxies (with Lfir, ^ 10 11 L0) and their 
Ultraluminous cousins (ULIRGs, Lfir > 10 12 L ). These sys- 
tems will also be accreted into the cluster environment along with 
their less-obscured and less-active population as the clusters grow. 
Indeed, mid-infrared (mid-IR) surveys of clusters have identified 
a populat ion of dusty star-bursts w hose activity increases with red- 
shift (e.g. lGeach et al.l2 006. 2008). However, these mid-IR surveys 
can miss the most obscured (and potentially most active) galax- 
ies which are optically thick in the rest-frame mid-IR. If they are 
present in clusters - even in relatively low numbers - such ac- 
tive galaxies will contribute significantly to the star formation rate 
(SFR) in these environments and the metal enrichment of the in- 
tracluster medium. Hence to obtain a complete census of the star 
formation within clusters we need to survey these systems at even 
longer wavelengths, in the rest-frame far-infrared (far-IR), corre- 
sponding to the observed sub-millimetre and millimetre waveband. 

Over the past decade or more there have been a number 
of studies of rich clusters of galaxies in the sub-millimetre and 
millimetre waveband s (e.g. ISmail et alj[l997l ; IZemcov et a i] |2007l : 
iKnudsen et alj|2008l) . Many of these studies were seeking to ex- 
ploit the clusters as gravitational telescopes to aid in the study 
of the distant sub-millimetre galaxy (SMG) population behind the 
clusters, while others focused on the detection of the Sunyaev- 
Zel'dovich (SZ) emission. Due to the limitations of current tech- 
nologies direct detection of millimetre sources is restricted to those 
with fluxes Snoo^m > 1 mJy, or equivalently galaxies with SFRs 
> 300MQyr _1 - much higher than expected for the general 
cluster populations based on optical surveys. Nevertheless these 
studies have serendipitously detected a number of cluster galax- 
ies, although these ar e either atypical (e.g. central cluster g alax- 
ies. lEdge et alj fl 999)) or are not confirmed members (e.g. iBestl 
l2002HWebb et alj2005l) . More critically, with few exceptions these 
studies have all focused on the central 2-3 arcmin of the clusters, 
where the SZ emission and lensing amplification both peak, and 
have not surveyed the wider environment of the cluster outskirts 
whe re much of the obsc ured star formation is likely to be occurring 
(e.g. lGeach etalj2006l) . The two excep tions are the wide- field sur- 
vey of the 2 ~ 0.25 cluster A2125 bv lWagg et all J2009b and the 
survey of a n overdense region of the COSMOS field bylScott et al.l 
(2008) and lAustermann et al.l d2009al) . lw"agg et alj J20091) detected 
29 millimetre galaxies across a ~ 25-arcmin diameter region cen- 
tered on A 2125 of which none are claimed to be members. How- 
ever, the only redshift estimates available are based on crude radio- 
to-millimetre spectral indices, which are sensitive to both dust tem- 



perature ^and^edsluftjBlah^ 

survey dScott et al.l2008T Austermann et al.l2009al) covered a num- 
ber of structures, including an X-ray detected z — 0.73 cluster and 
concluded that the statistical overdensity of sources was most likely 
due to the gravitational lensing of background SMGs by these fore- 
ground structures. 

To help to conclusively answer the issue of the obscured star- 
forming populations in distant clusters we have therefore under- 
taken a wide-field 1.1 -mm survey of the z = 0.54 rich cluster 
MS 0451.6-0305 (here after MS 0451-03) with the AzTEC cam- 
era dWilsonetai]|2008l) on the James Clerk Maxwell Telescope 
(JCMT). This panoramic millimetre survey can also take advan- 
tage of the significant archival data available for this well-studied 
region. In particular the panoramic imaging of MS 045 1 —03 from 
Spitzer and uniquely, the Hubble Space Telescope (HST), as well 
as extensive archival multiwavelength imaging and spectroscopy, 
which we employ for determining cluster membership of AzTEC- 
detected galaxies. Our survey probes the rest-frame 700 \im emis- 
sion of cluster galaxies - in search of examples of strongly star- 
forming but obscured galaxies - as well as identifying more lumi- 
nous and more distant examples of the SMG population. We can 
compare our millimetre search for cluster member s to the previous 
mid-IR survey of this cluster bv lGeach e t al.(2006) who uncovered 
a population of dusty star-forming galaxies which dominate the in- 
tegrated SFR of the cluster of 200 ± 100 M yr _1 . Our AzTEC 
map covers ~ 60 times the area of the 850 y im SCUBA observa - 
tions of the central region of MS 0451-03 dBorvs et alj |2004b), 
while the depth of a ~ 1.1 mjy is sufficient to identify ULIRGs 
individually and obtain stacked constraints on far-IR fainter cluster 
members. 

We describe our observations and the data reduction in ^2j 
present 1.1 -mm number counts, identify counterparts to the AzTEC 
galaxies, and use photometric techniques to determine cluster 
membership in i|3] ij4]contains our conclusions; individual AzTEC 
galaxies are presented in Appendix [A] Throughout this paper we 
use J2000 coordinates and ACDM cosmology with Qm = 0.3, 
Qa = 0.7 and Ho = 70 kms _1 Mpc _1 and all photometry is on 
the AB system unless otherwise stated. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 AzTEC Observations 

The AzTEC millimetre camera dWilson et alj|2008l) was installed 
on the JCMT during November 2005 and operated nearly contin- 
uously until February 2006. On the JCMT AzTEC has a field-of- 
view of ~ 18-arcmin 2 and a beam of 18-arcsec FWHM. For this 
project we observed a 0.23-deg 2 region centred on MS 0451—03 at 
04 h 54 m 10:l, 03°01'07'.'6, over a total of 24.5 hours between 2005 
December 14 and 2006 January 6. Observations were carried out 
in raster-scan mode with 1 200" long scans in e levation separated 
by 9" steps in azimuth (see Wil son et alj d2008l) for details of our 
standard raster scan strategy). Scans were made at a fixed speed of 
120"s _1 . Data from the turnarounds were excised for this analysis. 
Each observation took 35 minutes to complete and a total of 42 ob- 
servations were made of the field. Over the set of observations the 
zenith opacity at 225 GHz, T225, was monitored with the Caltech 
Submillimeter Observatory's opacity meter and ranged from 0.03 
to 0.17. The average opacity for the set of observations was 0.09, 
equivalent to an atmospheric column of precipitable water vapour 
of ~ 2 mm. 
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Figure 1. AzTEC S/N contours for the MS 045 1—03 field shown over the corresponding Subaru ft-band image; smoothed (dot-dashed) X-ray contours 
highlight the cluster centre. We also plot a 5-arcmin radius circle (dashed) illustrating the sub-region analysed in i]3.1l The solid black line encircles the 
region with exposure coverage ^ 50% of the maximum value, where sources are extracted. In this region AzTEC contours shown are at lcr intervals, starting 
at 2.5cr. 15 arcsec radius circles mark the 1.1-mm source detections labelled in order of decreasing S/N (Table|2j and squares highlight sources which are 
potential cluster members ( £|3.3.3t . MMJ 045413.35 is not shown because it lies outside of the area covered by the Subaru image. The 1.1-mm sources appear 
evenly distributed, with no obvious overdensity towards the cluster centre. 



The 42 individual AzTEC observations were reduced and 
combined using the publicly available AzTEC Data Reduction 
Pipeline V 1.0 - a set of IDL code that is optimised for detect- 
ing point sources in blank-field AzTEC surveys. This is the same 
pipeline that has be en used in p r evious blank-field A z TEC/ JCMT 
analyses including IScott et alj J2008t) . IPerera et all j2008h and 



lAustermann et al.l O009b ) 



Since the reduction of these observations follows the same 



procedure and use s the same code as previous AzTEC studies, we 
refer the reader to IScott et ai] d2008h for a thorough review of the 
reduction steps. Here we focus on the particulars of this field and 
the measured characteristics of the MS 045 1—03 map. 

Since we are interested in counterparts to our detected SMGs 
it is critical that we understand the pointing accuracy of the 
maps. In the reduction pipeline, fine corrections to the JCMT's 
pointing model are applied to each observation based on reg- 
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ular observations of a nearby QSO (a = 04 h 23 m 15°8, 5 = 
— 01°2 / 3 3'.'l) in the same m anner as described in IScott et all 
(2008) and Wilson et al . (2008). The residual astrometric error for 
the field is measured by stacking the AzTEC map at the positions 
of radio sources in the field (see i]2.2t and fitting for the centroid 
of the resulting image. The stacked radio sources are detected at 
~ 6cr and demonstrate an overall systematic astrometric shift in 
the AzTEC map of Aa = -3.0" ± 0.9", AS = 1.2" ± 1.5" with 
respect to the radio reference frame. This correction is negligible 
in Dec, but the significant RA offset is applied to the AzTEC data 
prior to further analysis and source extraction. 

The AzTEC data are flux calibrated as described in 
Wilso n et al. I d2008l) from nightly observations of Uranus over the 
JCMT run. The final calibration accuracy is ~ 10%. 

The primary products which come out of the AzTEC pipeline 
are a map representing the average filtered point source response, 
a map of the of the weight of each pixel in the sky flux map (rep- 
resenting the uncertainty in flux of each pixel), and a map of the 
signal to noise estimate in each pixel. All maps are made on the 
same 3" x 3" grid and all maps have been Wiener filtered to op- 
timize sensitivity to point sources. The central ~ 0.10 deg 2 of the 
AzTEC map is relatively uniform, with all pixels having ^ 50% of 
the peak weight. Sources are extracted in the central region of the 
maps ( i]3.1| > where the rms noise is ~ 1.1 mjy. Contours of 'signal- 
to-noise' (S/N), overlayed on the Subaru R-band image of the field 
are shown in Fig.[T] 

2.2 Multi-wavelength data 

Significant multiwavelength archival data exist for the 
MS 0451—03 field, which are summarised in Table Q] and 
described below. 

2.2.1 Spitzer MIPS imaging 

Geach et al. (2006) obtained Spitzer MIPS 24[i.m observations of a 
0.23-deg 2 area of MS 0451—03, excluding the central 25-arcmin 2 . 
The full det ails of the reduction and source extraction process are 
described in Ge ach et al. I d2006h . The central region of the cluster 
was part of Guaranteed Time Observations (Spitzer program 83) so 
these data were obtained from the archive and incorporated into the 
mosaic. The 5-a catalogue detection limit corresponds to 200 ixJy. 

2.2.2 Radio Imaging 

Archival observations of the MS 0451-03 field at 1.4 GHz 
were obtained using the National Radio Astronomy Observa- 
tory's (NRAC0) Very Large Array (VLA), combining 9 hours of 
data obtained in 2002 June in the VLA's BnA configuration with 
16 hours of A-configuration data taken in 2006 February (Project 
IDs AN109 and AB1199, respectively). The nearby calibrator, 
0503+020, was used to track amplitude and phase, with absolute 
flux and bandpass calibration set via 0137+331. The n ow-standard 
1.4-GHz wide-field imagin g approach was adopted l lOwen et"al] 
120051 : iBiggs & Ivisorj|2~006l) , using spectral-line mode '4' to ac- 
quire data with an integration time of 3 s ( 10 s in BnA) and using the 
IMAGR task in the Astronomical Image Processing System (MPS) 
to map out the primary beam using a mosaic of 37 images, with 17 

1 NRAO is operated by Associated Universities Inc., under a cooperative 
agreement with the National Science Foundation. 



more distant radio sources covered by additional facets. Several it- 
erations of self-calibration and imaging resulted in a noise level of 
11 /xjybeam -1 , witha2.3" x 1.8" synthesised beam at a position 
angle 0.0°. 

Sourc es and corresponding fl uxes are obtained from SEx- 
TRACTOR ( iBertin & Arnoutsl [19961) , using the S/N map for detec- 
tion and the flux map as the analysis image; sources are extracted 
where a minimum of 10 contiguous pixels (each is 0.16 arcsec 2 ) 
have S/N^ 2. The resulting catalogue has a 3-a flux limit of 51 
(jjy and is corrected for bandwidth smearing. We verify the statis- 
tical properties of our catalogue by comparing source counts with 
IBiggs & Ivisonl d2006t) . Notably, the source density is not signifi- 
cantly enhanced towards the cluster centre, allowing us to use the 
counts across the whole field in our statistical calculation of AzTEC 
counterparts ( i]3.2.2t . 



2.2.3 Optical and near-infrared imaging 

Observations of MS 0451—03 in the (/-band were taken with 
the MegaPrime camera on the Canada France Hawaii Telescope 
(CFHT) and reach a 3-a detection limit of 25. 1 mag in a 0.25-deg 2 
field. Standard reduction techniques were employed, customised 
for MegaPrime data (C.-J. Ma & H. Ebeling, p r ivate c ommunica- 
tion), details of which can be found in lDonovanl |2007). 

In addition we observed a 0.08-deg 2 area centered on the clus- 
ter core on 2007 October 09 through the z' filter with the Prime 
Focus Imaging Platform mounted on the 4.2-m William Herschel 
Telescope (WHT). A total integration time of one hour was ob- 
tained and reduced using standard techniques and calibrated using 
a Sloan Digital Sky Survey standard field, yielding a 3-er limiting 
magnitude of 24.2 mag. 

Finally, the BVRI imaging employed here comes from the 
PISCES survey (Panora mic Imaging and Sp ectroscopy of Cluster 
Evolution with Subaru; Kodama et al. 2005). Detection limits are 
listed in Table Q] In addition /f-band imaging to a 3-a depth of 
20.1 mag over 0.12-deg 2 centred on MS 045 1-03 was obtained 
from the Wide-Field Infrared Camera (WIRC; lwiTson et alj|2003h 
on the Palomar Hale telescope. Standard reduction methods were 
employed and are detailed in G. Smith et al. (in prep.). T hese 
data were used in th e previous analyses of lGeach et al. (2006) and 
Mora n et al. I d2007alTbh . 



2.2.4 Spitzer IRAC imaging 

The Spitzer InfraRed Array Camera (IRAC; Fazio et al. 2004) ob- 
servations of the 20' x 20' field centered on the MS0451 field were 
obtained as a part of the Cycle 5 General Observer (program 50610, 
PI: M. Yun) on 2009 March 18. Each tile was observed for a total of 
1500 seconds in each of the four IRAC bands (3.6, 4.5, 5.8, and 8.0 
Hm) in full array mode with a 15 position dither pattern with 100 
second exposures at each position. The new data are combined with 
the archival IRAC data (program 8 3) to produce the final mosaic 
images using CLUSTER GRINDER dGutermufh et al.ll2009T) . which 
is an IDL software package that utilizes standard Basic Calibrated 
Data (BCD) products from the Spitzer Science Centre's standard 
data pipeline. The angular resolution of the final mosaic ranges be- 
tween 2.0 to 2.5 arcsec depending on the observing band; the lim- 
iting depths are given in TableQ] 
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Table 1. Summary of the observations 



Filter 


Instrument 


Detection 


Reference 






Limit^ 




U 


CFHT - MegaPrime 


25.1 mag 


Donovan (2007) 


B 


Subaru - Suprime-Cam 


26.6 mag 


Kodama et al. (2005) 


V 


Subaru - Suprime-Cam 


25.8 mag 


Kodama et al. (2005) 


R 


Subaru - Suprime-Cam 


25.1 mag 


Kodama et al. (2005) 


I 


Subaru - Suprime-Cam 


24.2 mag 


Kodama et al. (2005) 


zf 


WHT - PFIP 


24 2 mag 


This work 


K 


Palomar Hale - WIRC 


20.1 mag 


G. Smith et al., in prep. 


3.6 |am 


Spitzer - IRAC 


23.0 mag 


This work 


4.5 \xm 


Spitzer - IRAC 


23.2 mag 


This work 


5.8 |am 


Spitzer - IRAC 


21.4 mag 


This work 


8.0 urn 


Spitzer - IRAC 


22.1 mag 


This work 


24 [am 


Spitzer - MIPS 


120 ^Jy 


Geach et al. (2006) 


1.1 mm 


JCMT - AzTEC 


3.6 mjy 


This work 


1.4 GHz 


VLA 


51 vJy 


This work 



a 1.1 mm is the 3.5cr minimum catalogue limit in the roughly constant noise region 
of the map; in the other bands we quote 3<r limits. 



2.2.5 Cataloguing 



To estimate redshifts photometrically we require that the photome- 
try in all filters samples the same emission from the source so an ac- 
curate spectral energy distribution (SED) can be built. To meet this 
requirement we degrade all the optical images to match the worst 
seeing - the ( V-band) in which the FWHM is 1 .54 arcsec. Each see- 
ing convolved image is astrometrically matched to the USNO cat- 
alogue with a set of unsaturated an d unblended stars spread across 
the frames. We use SEXTRACTOR dBertin & Arnouts| [l996*) on the 
i?-band image to detect objects with a minimum of ten adjacent 
0.2-arcsec pixels at least 1.5 a above the background to provide 
a source list and then use the APPHOT routine in IRAF to extract 
3-arcsec diameter aperture photometry at these positions in each 
convolved image. These measurements are then aperture corrected 
assuming a point source, to yield total magnitudes. We report 3-cr 
detection limits in Table Q] 



Survey limits at z=0.54 
Arp 220 scaled to 1.1mm 
M82 scaled to 1.1mm 



100 

Rest-Wavelength (^m) 



Figure 2. Catalogue detection limits (see Table [TJ in each filter for a 
z = 0.54 galaxy overlayed with Arp 220 and M82 SEDs ISilva et al J 1998h 
scaled to our 3.5<r 1.1 mm catalogue limit. We expect to easily detect coun- 
terparts at this redshift suggesting that unidentified SMGs are most likely to 
be high redshift background galaxies and not cluster members. 



The resolution of the IRAC images is significantly lower than 
the optical so for source extraction we consider these data sepa- 
rately, although an equivalent procedure is foll owed. Sources are 
detec ted on the 8 iim image with SEXTRACTOR teertin & Arnoutsl 
1996) and are required to have a minimum of 4 adjacent 0.9-arcsec 
pixels at least 2 a above the background. The 8 iim band is cho- 
sen so that we can use the IRAC colours to constrain counte rparts 
for otherwise unidentified SMGs ( ^3~231 and lVun et alj2008t) . This 
source list provides positions for the APPHOT routine in IRAF to 
extract 3.8-arcsec diameter fluxes. We also ensure that all SMGs 
which are identified in Submillimeter Array (SMA), radio or 24 
[im data and have IRAC counterparts in any bands are extracted, re- 
moving the requirement for an 8 txm detection. The extracted fluxes 
are corrected for the apertu re losses, employing the factors calcu- 
lated by the SWIRE team JSurace et aT1l2007t) . resulting in total 
magnitudes, which can be directly compared to our optical cata- 
logue. 



3 ANALYSIS AND RESULTS 

This study aims to identify 1.1 -mm detected ULIRGs in the z — 
0.54 galaxy cluster MS 045 1 —03. A priori, based on the space den- 
sity of ULIRGs at z ~ 0.5 Ee Floc'h et al .l2005h. and the overden- 
sity of LIRGs in clusters at this redshift iGeach et al.l2006|) . we ex- 
pect to find, at most, only a couple of cluster members in our survey. 
The large background population coupled with this small number 
of expected members makes this study challenging. Therefore, we 
focus on techniques to identify sources and confirm cluster mem- 
bership. Similarly to other SMG studies, a fraction of our sources 
are unidentified in the radio, mid-IR or optically. However, based 
on an Arp 220 SED, cluster members which are detectable above 
our 1.1 -mm flux limit are also expected to be brighter than the cata- 
logue limits in all other bands (Fig.[2]l. Therefore, any unidentified 
sources are likely to be part of the background population. 



© 2009 RAS, MNRAS OOO.fUflOl 



6 J. L. Wardlow et al. 



Table 2. AzTEC galaxies in the MS 0451—03 field in order of decreasing 
S/N. Source names correspond to AzTEC positions in the J2000 epoch. 
SMGs with robust multiwavelength counterparts ( £|3.2| and Table [5J are 
shown in bold. 





Source 


SNR 


Flux 


Noise 


Deboosted flux 








(mjy) (mjy) 


(mjy) 


1 


MMJ 045438.96-030739.8 


6.5 


8.3 


1.3 


64 -1.4 


2 


MMJ 045447.55-030018.8" 


5.8 


6.2 


1.1 


A — +1.1 

*■ -1.1 


3 


MMJ 045433.57-025204.0 


5.7 


7.7 


1.4 


~ .+1.4 
5 - 4 -1.5 


4 


MMJ 045431.56-025957.8 


5.1 


5.2 


1.0 


3 - 7 -l.l 


5 


MMJ 045421.55-030109.9 


5.1 


5.1 


1.0 


„ e + 1.0 


6 


* % % T t f \ 4 » 4 | m A 4~ \ 4~ \ ~V 4~ \ ~V 4~ \ S S 

MMJ 045417.49-030306.6 


4.9 


5.0 


1.0 


O e + 1.1 

d '°-1.2 


7 


MMJ 045413.35-031204.2'' 


4.6 


6.3 


1.4 


o ,1 + 1.4 
36 -2.0 


8 


A 1 mm I i \ A ^ A % ~\ ^7 ■> i\ "> 4\i\ A ~> 1-1 

MMJ 045412.72-030043.7 


4.5 


4.5 


1.0 


o n +l .0 


9 


MMJ 045345.31-030552.2 


4.4 


4.6 


1.0 


O C3 + 11 

2 ' 8 -1.4 


10 


MMJ 045407.14-030033.9 


4.3 


4.3 


1.0 


r, c + 10 

26 -1.4 


1 1 


MMJ 045358. 12-025233.3 


4.3 


4.3 


1.0 


n c + 10 

2 - 6 -i. 5 


12 


MMJ 045426.76-025806.5 


4.2 


4.2 


1.0 


2 ' 6 -1.4 


13 


MMJ 045356.09-02503 1 .4 


4.2 


5.4 


1.3 


2 - 6 -2.2 


14 


MMJ 045424.53-03033 1 .7 


4.1 


4.2 


1.0 


n ,1 + 1.0 

2 - 4 -1.5 


15 


MMJ 045328.86-030243.3 


4.1 


5.0 


1.2 


^-2.0 


16 


MMJ 045354.64-030004.0 


4.0 


4.0 


1.0 


r» i +1.0 

2 - 1 -1.6 


17 


MMJ 045431.35-025645.8 


4.0 


4.0 


1.0 


r, o + l 

2 - 3 -1.4 


18 


MMJ 045411.57-030307.5 


4.0 


4.1 


1.0 


23 -1.6 


19 


MMJ 045415.53-025 125.0 


3.9 


4.1 


1.0 


9 q+1.1 


20 


MMJ 045345.88-030440.0 


3.9 


4.0 


1.0 


2 - 2 -1.6 


21 


MMJ 045358.49-025601 .4 


3.9 


3.8 


1.0 


n o+l O 
2 - 2 -1.5 


22 


MMJ 045403. 10-025006.8 


3.9 


4.7 


1.2 


2 ' 2 -2.0 


23 


A /f~\ If T f~\ A iT C *1 A f f \ on^ 1 1 1 

MMJ 045357.46-030237. 1 


3.9 


4.0 


1.0 


n 1+1.0 


24 


MMJ 045420. 10-030658. 2 


3.8 


4.0 


1.0 


n o+l l 
2 - 2 -1.6 


25 


MMJ 045422. 17-030445.8 


3.8 


3.9 


1.0 


n n+10 
2 -0 1 7 


26 


MMJ 045349.69-025807.1 


3.8 


3.8 


1.0 


2.0+1° 


27 


MMJ 045421.17-030740.2 


3.8 


3.9 


1.0 


1 9 +1 1 


28 


MMJ 045345.06-025722.6 


3.7 


3.7 


1.0 


1 9+ 10 


29 


MMJ 045442.54-025455.0 


3.7 


4.5 


1.2 


1 6+ 14 

1 -°-1.6 


30 


MMJ 045444.78-030030.9 


3.7 


3.8 


1.0 




31 


MMJ 045411.73-025712.7 


3.7 


3.6 


1.0 


i-s±i:S 


32 


MMJ 045411. 17-031019.2 


3.6 


4.1 


1.1 


1.811;? 


33 


MMJ 045352.72-030130.9 


3.6 


3.7 


1.0 


2 0+ 10 


34 


MMJ 045454.20-025919.3 


3.6 


4.6 


1.3 




35 


MMJ 045359.27-030327.4 


3.5 


3.6 


1.0 


1 S+ 1 - 
L8 -1.6 


36 


MMJ 045340.09-030334.2 


3.5 


3.6 


1.0 


1 7 +1.0 
-"-■'-1.6 



a These SMGs were observed with the SMA and detected ( ^3.2. 1> . 

b These SMGs were observed with the SMA but not formally detected 

3.1 AzTEC catalogue 

Millimetre sources are identified as local maxima with S/N^ 3.5 
and are selected in the roughly uniform noise (a ~ 1.1 mjy, ~360- 
arcmin 2 ) region of the AzTEC signal-to-noise map (with weights 
of at least half of the maximum). In Table [2] we present these 36 
detections in order of decreasing S/N. Measured source fluxes and 
their errors are given by the peak pixel value and corresponding 
noise value respe ctively; deboosted fluxes are calculated follow- 
ing the method of IPerera et alj {2008). We determine the positions 
of sources on the sub-pixel level by calculating the centroid af- 
ter weighting pixels close to the brightest pixel according to the 
square of their fluxes. Fig. [T] shows the Subaru /?-band image with 
AzTEC contours of constant S/N and the detected 1.1 -mm sources 




1 2 3 4 5 6 7 

S (mjy) 



Figure 3. Differential number counts using deboosted fluxes from 
the AzTEC survey of MS 045 1 — 03 compared to those from the 
AzTEC/SHADES survey l Austermann et al. 2009tj). We see that across the 
whole area of our survey there is no strong overdensity of 1 . 1 mm sources in 
MS 0451-03 compared to the AzTEC/SHADES blank-field survey; hence 
it is apparent that the majority, but not necessarily the entirety, of the mil- 
limetre sources in the cluster field are background galaxies. However, in the 
central 5 arcmin radius of the map (~ 1.9 Mpc at z = 0.54) there appears 
to be a slight overdensity of sources. This could be due to a combination of 
the Sunyaev-Zel'dovich effect, gravitational lensing of background sources 
by MS 0451—03, and potentially a contribution from 1.1 mm cluster mem- 
bers. The solid line is the Schechter function fit to the number counts in the 
entire AzTEC map, and the dotted lines represent the survey limits - the 
lower and upper lines represent the whole AzTEC MS 0451—03 field and 
the central 5 arcmin radius, respectively. For clarity number counts in the 
central 5 arcmin radius of MS 045 1-03 and the AzTEC/SHADES survey 
are offset slightly in flux. 

ma rked. Based on the a nalyses of AzTEC data in lScott et alJ ( l2008T) 
and Per era et al. | d2008h we expect that 3^1 (8-10%) of the sources 
in our 3.5a catalogue are false detections with only ~ 0.5 and 
amongst those with S/N^ 4.5 and 5 respectively. 

To examine whether there is an overdensity of AzTEC galax- 
ies in MS 0451—03 compared to similar field surveys - poten- 
tially indicating a significant number of obscured ULIRG clus- 
ter members - we calculate the number counts in MS 0451—03 
compared to the AzTEC b lank-field survey of the SHADES fields 
dAustermann et alj |2009b ) . The source number counts and corre- 
sponding errors were ca lculated using the same bo otstrap sampling 
methods as described in lAustermann et ail d2009bl) . employing the 
AzTEC/SHADES best fit Schechter parameters as a prior. Due to 
the effects of flux boosting, our catalogue, which is limited at an 
apparent S/N ratio of 3.5, corresponding to a typical apparent flux 
of 3.6 mJy, actually contains sources fainter than this limit, allow- 
ing us to statistically constrain counts fainter than 3.6 mjy. We 
also calculate number counts in the central 5 arcmin radius of the 
AzTEC MS 045 1—03 map, by trimming the map to this area and re- 
peating the calculations with the sky model, deboosting, and com- 
pleteness estimates based on the full maps, since these are not ex- 
pected to change across the field. The number counts for the whole 
MS 0451—03 survey area, and the central 5 arcmin (which repre- 
sents a physical scale of ~ 1.9 Mpc at z — 0.54), compared to the 
AzTEC/SHADES survey are presented in Fig. [3] 
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Across the whole field MS 045 1 —03 does not exhibit a source 
excess at 1.1 mm compared to a blank field. This suggests that, as 
expected, there is not a dominant population of luminous obscured 
star-forming galaxies in MS 045 1—03. Nevertheless, it is possible 
that a small number of the sources are still cluster members. In- 
deed, our number counts analysis (Fig. [3} suggests there may be a 
small overdensity of sources in the central 5-arcmin radius region 
of the AzTEC map. Based on integral number counts and their er- 
rors at 1.1 mm down to 1 mjy, derived from bootstrap sampling 
the Posterior Flux Densit ies of sources (as described in detail in 
Austermann et al. (2009b)) a tentative overdensity is apparent at 
the ~ 1.6a level. If real this overdensity could arise from a com- 
bination of the Sunyaev-Zel'dovich effect, the gravitational lensing 
of background sources by the cluster, and potentially 1.1 mm clus- 
ter members. To discover whether any of this excess is caused by 
1.1 -mm bright cluster galaxies we must first accurately locate the 
counterparts to the 1.1 -mm emission. 

3.2 Identifying counterparts to AzTEC sources 

The large AzTEC beam (18-arcsec FWHM on the JCMT), coupled 
with the high spatial density of optical sources makes identification 
difficult unless precise positions are determined for the SMGs. The 
observed mm/sub-mm emission from SMGs represents rest-frame 
far-IR emission from dust-reprocessed starlight or AGN activity. 
Therefore, by expl oiting the far-IR-radio correlation (e.g. lCondonl 
ll992l ; lGarrettll2002l) and the low spatial density of radio sources, it 
is possible to employ deep (a ~ 10 p-Jy), high-resolution, interfer- 
ometric radio observations to deter mine accurate source positions 
for ~60-80 per ce nt of SMGs (e.g. llvison et al JI2002L [20051 120071 : 
Chapi n et alj2009h . Similarly, 24-p.m observations can be used for 
identification, a method which has proved useful in confirming ten- 
tative radio counterparts, or providing positions for a small number 
of the radio undetected SMGs (e.g. Ilvison et al.l 2004: Pop e et all 
l200d Ilvison et alj|2007l ; IChapin et alj|2009h . In addition, mid-IR 
colours have been used to isolate potential cluster SMGs - al- 
though this emission does not track far-IR outp ut and so it is an 
indirect indicato r of the luminous far-IR source dAshbv et al .120061 : 
IPope et alj2006l;lYun et alj2008h . 

An alternat i ve approach is submi l limetre interferome try (e.g. 
Ilono etal.ll2006l ; I Younger et ai] 120071 : ICowie et all 120091) . which 
is preferable to these traditional radio and mid-IR identification 
methods because it samples the same part of the Spectral Energy 
Distribution (SED) as our AzTEC imaging. However, sources 
generally require individual observations due to the small field- 
of-views of such instruments, commanding prohibitive exposure 
times to accurately locate all sources in a catalogue. 

We employ all three methods in this work, using millimetre 
interferometry where available, and radio and mid-IR imaging as 
an alternative. The details of each of these identification methods is 
discussed in the following sections. In total we identify 18 AzTEC 
counterparts (50% of the sample; Table [3}, of which 14 (39% of 
the total) also have optical counterparts (Table [4j. In Fig. [2] we 
demonstrate that at z ~ 0.5 our multi wavelength observations are 
deep enough to detect galaxies above our AzTEC flux limit and 
hence provide counterparts, assuming SEDs typical of Arp 220 or 
M82. Therefore, our unidentified sample is not expected to contain 
cluster members. Of the optically-unidentified SMGs, one lies out 
of the field of the Subaru, CFHT and IRAC observations, two are 
contaminated by nearby saturated stars and one is fainter than our 
detection limits at all wavelengths. We discuss the AzTEC galax- 



ies and corresponding counterparts on a source by source basis in 
Appendix lAl 

3.2.1 SMA detections 

The five brightest SMGs presented here were observed at 890 pm 
with the SMA in October and November 2007. On-source inte- 
gration times of 5-6 hours were employed, yielding maps with 
o"890(jm ~ 1.6 mjy. The SMA configuration resulted in a beam 
of ~ 3 arcsec. These data are discussed in detail in Wardlow et al. 
(in prep.). 

Of the five sources observed with the SMA MMJ 045438.96, 
MMJ 045447.55 and MMJ 045433.57 were detected with fluxes 
of 8.4, 5.5 and 8.3 mjy respectively. The remaining two targets, 
MMJ 045431.56 and MMJ 045413.35, were undetected, but have 
AzTEC 1.1 mm detections with SNR = 5.1 and 4.6 respec- 
tively. Unlike this MS 045 1—03 study, previous SMA observa- 
tions of AzTEC galaxies in the COSMOS field detected all seven 
targets \ Younger et alj|2007h . However, the faintest of the COS- 
MOS SMGs has deboosted 1.1 mm flux of 5.2 mjy, compared to 
3.7 and 3.6 mjy for MMJ 045431.56 and MMJ 045413.35 respec- 
tively. Therefore, based on the deboosted fluxes it is possible that 
MMJ 045431.56 and MMJ 045413.35 are too faint to be detected 
in our observations, although there are examples of bright SMGs 
that are undetected with the SMA (e.g. iMatsuda et~ai]|2007h . Al- 
ternatively, the lack of SMA counterparts could suggest that ei- 
ther the 890[J.m to 1.1 mm flux ratios of these sources are lower 
than expected (i.e. the dust is colder, or they are higher redshift 
than spectroscopically identified SMGs), that these galaxies ex- 
hibit extended far-IR emission on scales 3> 2 arcsec (2> 16 
kpc at z — 2; for example, from a merger-induced starburst), 
or that the AzTEC beam contains a blend of multiple SMGs. 
Our radio and SMA observations have similar resolutions, there- 
fore, extended or multiple radio counterparts can indicate the ex- 
tended or multiple nature of (sub)millimetre emission. Based on 
the radio emission it is likely that MMJ 045431.56 is composed 
of multiple (sub)millimetre sources, as discussed in Appendix lAl 
However, MMJ 045413.35 is more likely to be a single resolved 
(sub)millimetre source, with a lower 890-to-1100 pm ratio than 
expected. These issues are discussed further in Wardlow et al. (in 
prep.). We use the positions of the three SMA detected SMGs in 
the following analysis. 

We note that the proposed radio identifications of 
MMJ 045438.96 and MMJ 045433.57 are coincident with their 
SMA positions. The other SMA detected galaxy, MMJ 045447.55, 
lies ~ 9 arcsec from an elliptical galaxy with bright radio and 
mid-IR emission, which, due to the high fluxes is formally a 
'robust' identification, but is unlikely to be the true source of the 
millimetre emission. MMJ 045431.56 and MMJ 045413.35 which 
were targeted but not detected with the SMA are identified through 
radio and mid-IR counterparts respectively. 
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Table 3. Radio and 24 |am counterparts of AzTEC galaxies in MS 0451—03. All matches within 10-arcsec are listed, those secure identifications at 1.4 GHz or 24 |am with P ^ 0.05 are shown in bold and tentative 
associations (0.05 < P SC 0.10) are also presented. Images and a discussion of each source are presented in Appendix |A| 

h 

Source 1 .4GHz position^ Si.4GH2^1 l.lmm-1.4GHz Pi.4GHz 2A\aa position^ S , 24mn£ l.lmm-24|jm P24 ]xm 24mm-1.4GHz^ ^ 

RA Dec Separation RA Dec Separation Separation O 

(J2000) (uly) (arcsec) (J2000) (ijjy) (arcsec) (arcsec) g. 



1 


MMJ 045438.96/ 


04 h 54 m 39!01 


-03°07'38'/5 


98.0 


1.51 


0.003 


04 h 54 m 38!97 


-03°07'38"0 


157 


1.3 


0.023 


0.79 


2 


MMJ045447.55 - £ i 


04 h 54 m 47?46 


-03°00'19'.'2 


< 51 


- 


- 


- 


- 


< 120 


- 


- 


- 




04 h 54 m 48?14 


-03°00'14'.'6 


101 


9.78 


0.032 


04 h 54 m 48?12 


-03°00'14'.'3 


3396 


9.56 


0.009 


0.46 


3 


MMJ 045433.57^ 


04 h 54 m 33!55 


-02°52'04"6 


93.5 


0.71 


0.0008 


04 h 54 m 33!51 


-02°52'06"8 


392 


3.00 


0.020 


2.28 


4 


MMJ 045431.56^, 


04 h 54 m 31!69 


-03°00'07"1 


77.8 


9.50 


0.035 


04 h 54 m 3f?13 


-02°59'51'.'6 


485 


8.89 


0.076 


17.62 


5 


MMJ 045421.55 


04 h 54 m 21!66 


-O3°01'08"5 


84.3 


2.21 


0.006 


04 h 54 m 21!40 


-O3°01'07"5 


209 


3.29 


0.046 


4.09 


6 


MMJ 045417.49 


04 h 54 m 17 s 27 


-O3°03'03"2 


189 


4.74 


0.010 


04 h 54 m 17?24 


-03°03'03'/7 


217 


4.76 


0.074 


0.58 


7 


MMJ 045413.35 


04 h 54 m 13!36 


-03°11'58"8 


78.4 


5.44 


0.021 


04 fh 54 m 13!38 


-03°11'58'.'9 


333 


5.31 


0.056 


0.29 


$ 


MM I 04^41? 7? 


rt4 h ^4 m 1'? s 79 


— uo uu ^j. y 


58.3 


1.03 


0.002 






s 1 9ft 








9 


MMJ 045345.31 


04 h 53 m 45!29 


-O3°05'53"5 


92.3 


1.30 


0.002 


04 h 53 m 45!36 


-O3°05'52"8 


1393 


0.86 


0.0004 


1.24 






04 h 53 m 45?07 


-03°05'56{'6 


76.9 


5.62 


0.022 












5.71 


10 


MMJ 045407.14 


04 h 54 m 07!08 


-03°00'37"2 


60.9 


3.37 


0.013 


04 h 54 m 07!07 


-03°00'37"5 


1120 


3.76 


0.007 


0.39 


12 


MMJ 045426.762. 


04 h 54 m 26?86 


-02°58 / 08 / /0 


< 51 










< 120 








15 


MMJ 045328.862 






~ 51 






04 h 53 m 28!83 


-03°02'45"9 


150 


2.54 


0.040 




17 


MMJ045431.35 _ 


04 h 54 m 31!71 


-02°56'45"0 


89.3 


5.41 


0.019 






< 120 








18 


MMJ 04541 1.57^ 






~ 51 






04 h 54 m ll!66 


-O3°03'08"0 


185 


1.55 


0.016 




26 


MMJ 045349.69 






< 51 






04 h 53 m 49!80 


-02°58'13"9 


530 


7.03 


0.049 




27 


MMJ045421.17J 






~ 51 






04 h 54 m 21!14 


-03°07'46"1 


464 


5.92 


0.045 




28 


MMJ 045345. 06 J 






~ 51 






04 h 53 m 45?33 


-02°57'18'.'l 


226 


6.00 


0.095 




29 


MMJ 045442.54* 


04 h 54 m 42?62 


-02°54'48 / .'2 


< 51 










< 120 









a 1.4 GHz positions have typical uncertainties of ~ 0.6-arcsec. 
^ 1.4 GHz fluxes have typical errors of 17 \xjy. 
24 |jm positions have typical uncertainties of ~ 1.8-arcsec. 
^ 24 torn fluxes have typical errors of 40 |xly. 

e 24nm-1.4GHz separations consistent with the same counterpart are italicised. 
' The radio and mid-IR positions agree with the SMA position ( i|3.2.U 

9 MMJ 045447.55 is identified from SMA observations (Wardlow et al., 2009b), but a bright elliptical galaxy ~ 9.6 arcsec from the AzTEC position is also a formal radio and mid-IR identification. We consider 

the radio and mid-IR galaxy to be a chance association - at P = 0.05 we expect 1-2 chance associations in our catalogue. 
® ^ There is an additonal 176-|aJy 24-|am source, with P = 0.17, coincident with the radio counterpart of MMJ 04543 1.56. It is considered related to the AzTEC detection due to the radio identification, 
o ' MMJ 045426.76 and MMJ 045442.54 have no radio or mid-IR counterparts but are identified on the basis of their IRAC colours lYun et alj2008l) . 
70 3 There is a ~ 3cr radio peak coincident with the 24- [jm counterpart. 
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3.2.2 Radio and mid-IR identifications 

In this work we use both 1.4-GHz VLA and 24- nm Spitzer MIPS 
observations to locate the AzTEC galaxies in our sample. We iden- 
tify SMG counterparts independently at radio and mid-IR wave- 
lengths before comparing these for each AzTEC source. In our 
maps the surface density of radio galaxies is lower than the mid-IR, 
the positional accuracy is greater (~ 0.6 arcsec compared to ~ 1.8 
arcsec), and the link from the radio emission to the far-IR is tighter. 
Therefore, if the 24 [im and 1.4 GHz identifications disagree we 
consider the radio position as the more reliable counterpart. 

To ensure that no genuine associations are missed we search 
up to 10 arcsec from each AzTEC position. This corresponds 
to a ~ 3<t search radius for the lowest signal-to-noise ratio 
(SNR) sources, where a is the error on AzTEC position, given by 
O.GFWHMfSNR) " 1 , and the FWHM is of the instrument beam 
dlvison et alj|2007l) . We reject counterparts with more than 10 per 
cent probability of be ing chan ce associations using the P-statistic o f 
iDownes et ai1il986h (see also lPope et alj200o1 ; llvison et al.l2007h . 
Counterparts found with P < 0.05 are considered robust and those 
with 0.05 < P < 0.1 tentative. We present all identifications in 
Table H] 



3.2.3 IRAC identifications 

lYun et al.l (|2008) studied SMGs securely identified with SMA, ra- 
dio, or 24 urn data and found that they typically have redder IRAC 
colours than the submillimetre-faint foreground galaxy population. 
They proposed a selection criteria for SMG counterparts, based on 
IRAC colours and found that the SMA detected galaxies in their 
sample, whether radio and mid-IR identified or not, were all re- 
covered by this method. We verify that this selection criteria for 
our SMA, radio and 24 ixm-identified SMG counterparts and use 
it to search for counterparts to our otherwise unidentifi ed AzTEC 
source s and identify two further galaxies. Unlike the lYun et all 
(2008) SMA-detected SMGs, MMJ 045447.55, our optically, radio 
and 24(j.m faint, but SMA-detected SMG, is not detected in any of 
the IRAC wavebands to the depth of the survey listed in TableQ] 



3.3.1 Simple photometric redshift analysis 

In Fig. g] we plot 524nm/Si.imm versus SiAGHz/Si.imm for 
AzTEC MS 045 1-03 SMGs. F or comparison we plot 850 \im 
sources from llvison et ail d2007h . extrapolated from the observed 
850 [im to the equivalent 1.1mm flux assuming a power law spec- 
trum of the form Su ~ u 3 ' 5 . As expected both sets of SMGs lie 
broadly within the same region of colour-colour space, adding con- 
fidence to our detections and identifications. The redshift tracks 
of the local star-forming galaxy Arp 220, and the higher r edshift 
HR10 (z = 1.44) (based on the SEDs of Isilva et al] ( 119981) ) sug- 
gest that the AzTEC galaxies have redshi fts of 1 < z < 3.5 , 
in concordance with the SMG population ([Chapman et al J 120050 . 
However, potential degeneracy between redshift and dust tempera- 
ture in the templates makes this method unreliable for identifying 
cluster galaxies. 

We next determine whether any of the SMGs are likely to lie 
at low redshift and hence p otentially be memb ers of MS 0451—03 
using the BzK selection of Dad di et al. to separate z > 1.4 

galaxies from those at z < 1. Of the ~ 70 SMGs with spectro- 
scopic redshifts presented in Chapm an et alj d2005h ~ 80% have 
z > 1.4, suggesting that the BzK selection should enable us 
to remove the majority of background sources from our sample. 
Fig.[5]shows the optically identified SMGs in MS 045 1-03 which 
are co vered by our B, z and K imaging, in addition to Dad di et all 
(2004) selection criteria for high redshift (2 > 1.4) galaxies; SMGs 
with z p hot < 1.4 are highlighted (^333}. It is clear from Fig. f5]that 
from a BzK selection alone, at least two SMGs are low redshift 
and therefore potential cluster members; several additional galax- 
ies lie close to the border or have photometric limits which could 
place them in the low redshift region. Therefore, we next carry out 
a full photometric redshift an alysis of the whole sample using Hy- 
PERZ jBolzonella et alj|200fjh to calculate redshifts for the SMGs 



3.3 Redshifts of MS 0451-03 AzTEC galaxies 

Our study seeks to identify potential millimetre sources in the clus- 
ter population of MS 0451—03. There are spectroscopic observa- 
tions of 1639 galaxies within our field, but nevertheless none of 
the AzTEC galaxies have been spectroscopically observed. There- 
fore, we use photometric methods to separate any potential AzTEC 
detected cluster members from 'typical' z ~ 2 SMGs. We first 
apply two simple colour tests - the advantage of these is that it 
is easy to understand the biases in the sample - before apply- 
ing a more sophisticated photometric redshift analysis. We use 
the 5 , 24nm/S'i.i mm and Si.4GHz/Si.imm flux ratios to estimate 
the redshifts of identified AzT EC galaxies and th en also consider 
the BzK selection criteria of Daddi et al. (2004) to separate the 
higher redshift (z > 1.4) SMGs from potential cluster members. 
Finally, we test whether these conclusions are reliable by using 
the spectrosc opic redshif t s and multi-wavelength photometry of 
the SMGs in lBorvs et al. (2005), and apply these findings to our 
AzTEC MS 0451—03 sources to obtain photometric redshifts for 
our galaxies. 
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Table 4. Optical and near-IR photometry for the detected SMG counterparts with derived photometric redshifts. Potential cluster members are shown in bold (J533J- 



O 



1 

3 

4 

5" 

8 

10 

12 

15 

17 

18 

26 

27 

28 

29 



Source 

MMJ 045438.96 
MMJ 045433.57 
MMJ 04543 1.56 
MMJ 045421.55 
MMJ 0454 12.72 
MMJ 045407. 14 
MMJ 045426.76 
MMJ 045328.86 
MMJ 045431.35 
MMJ 04541 1.57 
MMJ 045349.69 
MMJ 04542 1.1 7 
MMJ 045354.06 
MMJ 045442.54 



U 

> 25.1 
26.92±0.17 
24.04±0.04 
24.92±0.08 

> 25.1 
22.83±0.03 

> 25.1 

> 25.1 
24.45±0.05 
24.83±0.05 
26.41±0.14 
25.42±0.10 
25.61±0.08 
25.00±0.09 



B 

25.89±0.04 
25.60±0.04 
24.20±0.02 
24.50±0.03 

> 26.6 
22.23±0.02 
27.32±0.1 1 

> 26.6 
24.37±0.03 
24.65±0.03 
26.25±0.06 
25.51±0.04 
25.19±0.03 
26.16±0.05 



V 

24.03±0.04 
24.96±0.05 
23.45±0.03 
23.37±0.03 

> 25.8 
21.69±0.02 
27.71±0.57 

> 25.8 
24.23±0.03 
24.23±0.03 
27.72±0.38 
24.46±0.04 
24.83±0.05 
25.71±0.12 



R 

23.29±0.03 
24.02±0.03 
22.71±0.02 
22.31±0.02 

> 25.1 
21.23±0.02 
25.53±0.08 

> 25.1 
23.72±0.03 
23.53±0.03 
26.04±0.10 
23.56±0.03 
24.68±0.04 
24.61±0.04 



/ 

22.77±0.03 
23.35±0.04 
21.97±0.03 
21.44±0.03 

> 24.2 
20.80±0.02 
23.83±0.05 

> 24.2 
23.31±0.04 
22.82±0.03 
24.58±0.09 
23. 11 ±0.04 
24.35±0.07 
23.86±0.05 



23.75±1.10 

21.22±0.05 
21.31±0.09 

> 24.2 
20.64±0.04 
26.80±4.24 

> 24.2 
25.94±1.40 
22.68±0.11 
24.35±0.40 
23.68±0.22 
24.36±0.34 

> 24.2 



K 

19.88±0.11 
19.27±0.07 
18.32±0.06 
18.81±0.12 

> 20.1 
19.67±0.08 
22.39±0.34 

> 20.1 
21.13±0.17 
22.27±0.37 

> 20.1 
18.82±0.07 
21.63±0.43 
20.65±0.13 



3.6|jjn 

20.78±0.05 
20.41±0.03 
20.18±0.03 
19.70±0.02 
19.48±0.01 
20.78±0.05 
21.95±0.13 
21.08±0.06 
22.21±0.17 
22.37±0.20 



4.5 |jm 

20.37±0.03 
20.05±0.02 
20.14±0.02 
19.96±0.02 
19.03±0.02 
20.37±0.03 
21.77±0.10 
20.83±0.04 
21.97±0.12 
21.78±0.10 



5.8|am 

19.92±0.09 
19.72±0.08 
20.15±0.11 
19.43±0.06 
18.61 ±0.03 
19.92±0.09 
21.37±0.36 
20.48±0.15 
22.37±1.21 
21.41±0.37 



8|am 

20.44±0.07 
20. 11 ±0.05 
20.85±0.10 
20.54±0.08 
19.40±0.03 
20.44±0.07 
21.25±0.15 
21.48±0.19 
22.07±0.33 
21.80±0.25 



z p hot 



3.42 
2.55 



20.07±0.02 19.83±0.02 19.15±0.05 20.52±0.08 
21.59±0.10 21.42±0.07 20.49±0.16 21.31±0.16 
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a This is the combined photometry for C 1 and C2 as discussed in £|3.3.3I 
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Figure 4. S^m/Si. 1mm VerSUS S'l,4GHz/5'l.lmm for SMGs in the 
MS 0451—03 field. We differentiate between sources with robust and ten- 
tative detections (0.05 < P ^ 0.10), and high light SMGs which satisfy 
the BzK selection criteria of Daddi et al. 12004) (Fig. [5). Sources without 
both mid-IR and radio detections are plotted at the 3cr detection limit of 
the respective catalogues (Table Q}. For comparison we also plot SHADES 
identifications ilvison et al. 2007), converted to expected observed 1.1 mm 
fluxes as discussed in the text; error bars of this sample are omitted for 
clarity. Both sets of SMGs lie broadly within the same region of colour- 
colour space, providing confidence in our identifications. Redshift tracks 
of th e 2 = 1.4 SMG H R 10 and local starburst Arp 220 (based on SEDs 
from ISilva et all d 19981) ) suggest that the AzTEC sourc es generally have 
1 < z < 3.5, in concordance with[Chapman et al. 120(3). 



3.3.2 Photometric redshifts of SMGs 

Having located the MS 045 1—03 SMGs we now need to test 
whether any are potential cluster members. Photometric redshifts 
have been calculated for several sets of SMGs using various codes 
and spectral templates , generally designed f or use on 'normal' low 
redshift galaxies (e.g. IClements et al.l l2008). However, SMGs are 
located at high redshifts and powered by dusty starbursts which 
can be contaminated by AGN meaning that the stellar templates 
derived from low redshift galaxies may not be appropriate. 

Typically studies of SMGs do not have both spectroscopic 
and photometric information so testing of photometric redshift 
results is difficult. Of ten results from codes such as HYPERZ 
jBolzonella et alj|2000l) or IMPZ jBabbedge et alj|2004l) are com- 
pared to those from other, apparently cruder, methods of redshift 
esti mation - such as the radio to far-IR -(sub)-mm spectral indices 
(e.g.lCarilii&Yunl 19991 ; lYun & Carillil2002l; lAretxaga et alj2003l : 
Clem ents et al.ll2008t> . There are two main problems with this ap- 
proach - firstly, the comparison is made with a set of results which 
themselves are not known to be correct, and secondly, the errors 
on the comparison redshifts are typically large, such that a gen- 
eral agreement can be confirmed but nothing more. In cases where 
spectroscopic information has been o btained for a su b-sample of 
galaxies typically they are small (e.g. iPope et al.ll2005h . 

In this work we expand on previous analyses to examine how 
well redshifts can be constrained for a spectroscopically confirmed 
sample of typical SMGs (median z ~ 2.2). We employ the sam- 



Figure 5. B — z vs. z — K colour-colour plot of the optically identified 
AzTEC galaxies in MS 0451—03 which lie within the field-of-view of our 
B, z and K observations. The lines of separation between passive and star- 
forming 2 > 1.4 BzK galaxies (pBzK and sBzK respectively), z < 1.4 
galaxies, and stars I Daddi et al. 2004) are shown. Two SMGs occupy the 
2 < 1.4 region, suggesting they are potential cluster members. We dis- 
tinguish between galaxies with z p f lot > 1.4 and z v ^ ot < 1.4 based on 
photometric redshifts calculated in i|3.3.3| The BzK selection and full pho- 
tometric analysis broadly agree as to the high and low redshift samples giv- 
ing more confidence in our ability to isolate millimetre sources in the clus- 
ter from the dominant background population. As expected, by this criteria 
most SMGs are deemed to be actively star-forming. 



pie of 12 SM Gs from Borvs et al. 1 20051) w hich have spectroscopi c 
redshifts from lChapman et al] d2003Ll2005l) and lPope et all d2008l). 
Photom etry in 12 bands from U to 8 \im is presented in lBorvs et al] 
d2005l) . We note that SMM J 123622.65 has a rest-frame UV/optical 
spect rum with z = 2.466 (Chapman et al. 2005; Swinban k et al] 
2004) but mid-IR spectroscopic analysis suggests z — 1.79 ± 0.04 
(Pop e et al . 2008), therefore, we exclude it from the statistical anal- 
yses. We also exclude SMM J 1 237 1 2.05 from the sample as it is un- 
detected at wavelengths shorter than 2 [im and the power-law shape 
at longer wavelengths means that th e redshift cannot be der ived. 

We use the HYPER^] package dBolzonella et al.l200oh for our 
photometric redshift estimates. The program calculates expected 
magnitudes in the observed filters for given SEDs of a library of 
model star-formation histories, with different ages, reddening and 
redshifts; the observe d and expected magnitude s are compared in 
each filter. We use the Bruz ual & Charload 19931) spectral templates 
provided with the HYPERZ package which represent star-formation 
histories resulting in SEDs which match local ellipticals (E), Sb, a 
single burst (Burst) and a constant Star-formation rate (Im). Red- 
shifts between and 7 are considered. SMGs are known to be dusty 
systems, therefore, we i nitially allow r edden ing of ^ Ay ^ 5, 
in steps of 0.2 using the lCalzetti etail j2000h reddening law. Ages 
of the galaxies are required to be less than the age of the Universe 
at the appropriate redshift. 

In Fig. [6] we compare the spectroscopic with our photomet- 



We use HYPERZ version 10.0 (http://www.ast.obs-mip.fr/users/roser/hyperz/). 
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Figure 6. Spectro scopic versus photometric redshifts for the SMGs in 
iBorvs etalJ 120051) . Gala xy SMMJ 123622.65 has rest-frame UV/op tical 
spectra with z = 2.466 IChapman et alj|2005t ISwinbank et ai]|2004) but 
mid-IR analysis places it at z = 1.79 ± 0.04 IPope et alj2008l) . therefore, 
SMMJ 123622.65 is plotted at each redshift. Utilising 12-band photometry 
results in constraints that are sufficient for separating high-redshift SMGs 
from those at z ~ 0.5. 



ric redshifts for the iBorvs et alj J2005t) galaxies; error bars shown 
are the HYPERZ quoted 99% confidence intervals which represent 
the la limits more realistically than the HYPERZ quoted 68% in- 
tervals for this sample. The average scatter in Az = z p hot — z spcc 
is 0.27 ± 0.21, although this is smaller for the z < 1.5 SMGs 
(Az — 0.15 ± 0.13) than for the more distant sample (Az = 
0.32 ± 0.22), demonstrating that when testing photometric red- 
shifts of these galaxies it is important to consider a comparison 
sample with the same redshift distribution as the population under 
study. Our results are unchanged if we restrict Av to less than 2.5 
although a tighter limit than this does affect the results - increasing 
Az. 

AGN contamination in the 8 \im IRAC filter is possible for 
SMGs (Hainline, L. J. et al. in prep.), therefore we repeat the pho- 
tometric redshift calculations as above but excluding the 8 \im data 
from the analysis. For all the SMGs Az n o&\i.m ^ Az a nbands yield- 
ing an average scatter Az — 0.23 ± 0.22 for the high redshift 
(z > 1.5) galaxies - consistent with the view that the 8 \im fluxes 
are contaminated. 

We also find no change if we allow a full range of templates: 
E, SO, Sa, Sb, Sc, Sd, Burst, and Im spectral types, and minimal 
change if only we restrict the templates to only Burst and Im models 
- despite the fact that when the full range of templates is allowed 
most of the galaxies are best-fit by the burst and elliptical models. 
We therefore find it unlikely that the addition of any further pure- 
stellar templates will improve the redshift accuracy. 

In this work we are interested in separating low redshift (z ~ 
0.5; potential cluster member) SMGs from the typical high redshift 
(z ~ 2) population, and it is clear from Fig. [6] that this is possible 
based on the photometric redshift resolution. 

3.3.3 Photometric redshifts of the AzTEC galaxy sample 

We use HYPERZ to calculate the redshifts of the SMGs in the field 
of MS 0451—03 from our optica l, near-IR and mid-IR photome- 
try. The lBruzual & Chariot! J 19931) spectral templates included with 
HYPERZ are used, and, as discussed in S|3.3.2| we consider only the 
Burst, Elliptical, Sb and constant star-formation rate models, limit 
reddening to ^ Av ^ 5, and redshifts to ^ z ^ 7. Bright, re- 



solved galaxies with uncharacteristically high primary photometric 
redshifts are considered to lie at the calculated secondary solutions. 

The redshift estimate of each counterpart is presented in Ta- 
ble [4] and in Fig. [7] we show the SMG photometric redshifts ver- 
sus AzTEC 1.1 -mm fluxes. Although, most of the AzTEC sources 
are high redshift background galaxies, two - MMJ 045421.55 and 
MMJ 04543 1.35 - are possible cluster members, and are discussed 
further below. W e note that as reported by other authors (e.g. 
IPope et al.ll2005l) there appears to be a weak trend of redshift with 
millimetre flux - the brightest millimetre galaxies lie at higher 
redshifts. Our sample has a median photometric redshift of 1.2 
and an interquartile r ange of z = 0.7-1. 9 - lower than spectro- 
scopic studies (e.g. IChapman et al]|2005l (z) = 2.2) but simi- 
lar to other photometric studies (e.g. (z) = 1.4; Iciements et all 
2008). We note that if we exclude the potential cluster members 
MMJ 045421.55 and MMJ 045431.35 from this analysis the me- 
dian redshift of the field SMGs in this study is 1 .8 and the interquar- 
tile range is z = 0.9-1.9. 

Although in &|3.3.2| we showed that the exclusion of 8 p.m in- 
formation can sometimes improve the accuracy of photometric red- 
shift estimates of SMGs, the redshifts reported in Table [4] include 
the 8 |J.m photometry. This is because for the galaxies with only 
IRAC detections the resulting lack of information about the loca- 
tion of the 1.6iim stellar peak means that only weak redshift con- 
straints are possible if we remove the 8 \im photometry. Critically, 
the inclusion or exclusion of the 8 |J.m information does not affect 
MMJ 045421.55 and MMJ 045431.35 which still both have photo- 
metric redshifts consistent with the cluster. 

In Fig.|5]we showed the B — z versus z — K colour-colour plot 
for AzT EC galaxies in M S 0451-03 and used the BzK selection 
criteria lDaddi etal.ll2004l) to separate galaxies above and below 
z = 1.4. We similarly group the SMGs based on the full photo- 
metric analysis from HYPERZ and find that both methods broadly 
agree. In Fig. [7] we show the redshifts and AT-band magnitudes of 
MS 045 1—03 SMGs in comparison to t he spectroscopic sample of 
SMGs examined by ISmail et al] (120041), an d local ULIRGs from 
iKim et al.l fc002l) and lStanford et alj d2000). The apparent A"-band 
magnitudes of our AzTEC sample are consistent with their esti- 
mated redshifts, when compared with spectroscopic SMG surveys 
and local ULIRGs, suggesting that our photometric redshifts are 
reasonable. The general agreement between our photometric red- 
shift analysis, BzK and A^-band magnitudes supports the derived 
redshifts of the AzTEC sources. 

The potential cluster members, MMJ 045421.55 and 
MMJ 04543 1.35, are discussed in detail here. The other identified 
SMGs are examined in Appendix lAl 

MMJ 045421.55 MMJ 045421.55 is identified through radio 
emission 2.2 arcsec from the AzTEC centroid, which lies between 
two optical galaxies (2.1 and 1.6 arcsec from the northern (CI) 
and southern (C2) galaxies respectively). There is also an IRAC 
source at the location of the radio emission, which appears slightly 
extended towards C2 (Fig. lAlb . The three possible explanations 
for this system are: CI and C2 are interacting and both millimetre 
bright; C2 is the millimetre counterpart; or neither CI nor C2 are 
responsible for the millimetre emission. Each of these possibilities 
is discussed below. 

Individual photometric analysis of CI and C2 (excluding the 
IRAC information) suggests that they are both potential cluster 
members, and with a separation of 2.6 arcsec - corresponding 
to ~ 17 kpc at z = 0.54 - they could be in the early stages 
of a merger. Although high-resolution (0.1 arcsec) archival HST 
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Figure 7. The left-hand panel shows photometric redshifts of SMGs against the 1.1 mm AzTEC fluxes. Two SMGs - MMJ 04542 1.55 and MMJ 045431.35 
- have photometric redshifts consistent with being cluster members and are highlighted; we discuss these two sources in deta il in i]3.3.3l In the right-hand 
panel we sho w redshift versus if -b and magnitudes of M S 0451—03 SMGs compared to the spectroscopic sample presented in Smail et al. 12004), and local 
ULIRGs from lKimet a?] J2002h and lStanford et al l feOOOl) . Both datasets occupy the same parameter space suggesting it is unlikely that any of our photometric 
redshifts are extreme outliers. 



F814W imaging shows no evidence of disturbance (Fig. IAU , the 
expected tidal tails can be low surface brightness features, mak- 
ing them difficult to detect. Such interactions are widely known to 
trigger dusty starbursts in which the radio emission appears to be 
located between the optical nuclei (e.g. systems similar to VV 114; 
iFraver et alJl999tlLe Floc'h et alj2002tllono et al]2004h . 

With the aim of measuring redshifts we targeted CI and C2 
with the ISIS long-slit spectrograph on the WHT during service 
time in 2009 February. The total integration time was one hour in 
a — 1.1-arcsec seeing and standard reduction techniques were 
employed. A faint continuum was observed (the two targets are 
blended), but no features suitable for redshift measurement were 
detectable. A possible faint (~ 2a) emission feature at 7852 A 
is visible, which, if real, is most likely to be [OIII] A5007 A at 
2 = 0.568 - placing the galaxy in a small group just behind the 
cluster. However, the feature is tentative and inconclusive. 

On the assumption that CI and C2 lie at the same redshift, 
and that the millimetre, radio and IRAC fluxes are emitted from a 
merging system as described above, we combine the optical fluxes 
from CI and C2. The optical to near-IR photometric redshift of the 
whole system is z — 0.501q'iq which, as shown in Fig. [7] agrees 
with the z versus K ULIRG trend. Therefore, if C 1 and C2 are an 
interacting system, MMJ 045421.55 a possible cluster member. 

The radio counterpart lies only 2.6a from C2 (compared to 
3.5a from CI), and the IRAC emission of MMJ 045421.55 appears 
extended towards C2. Therefore, in the situation where CI and C2 
are unassociated we find it most likely that the C2 is the counter- 
part to the IRAC, radio and millimetre flux. In this case we ob- 
tain a photometric redshift of z = 0.511 o 'q5 - once again placing 
MMJ 045421.55 in the region of MS 0451-03. 

Based on the intrinsic faintness of many SMG counterparts it 
is feasible that the identified radio and IRAC emission arises from 
an optically faint background galaxy, unrelated to CI and C2, and 
undetectable in our observations. Previous SMG surveys have con- 
fused low redshift galaxies with the source of sub-mm emission due 



to the lensing of the background source (e.g. lChapman et al.l2002l) . 
In these situations it is only possible to distinguish between a lensed 
background galaxy and a foreground cluster member through the 
detection of the faint optical counterpart in deeper imaging, or by 
the detection of CO emission lines. 

MMJ 045431.35 The counterpart to MMJ 04543 1.35 is securely 
identified 5.4 arcsec from the AzTEC centroid through its radio 
emission. The source is detected in all our Spitzer and ground- 
based imaging and resolved in the HST F&14W image into merging 
galaxies with centroids separated by ~ 1.6 arcsec (~ 10 kpc at 
z = 0.54) and tidal tails between them (Fig. IAU . We suggest 
that a dusty starburst triggered by the interaction between the two 
galaxies is causing the millimetre emission from this system. The 
photometric redshift is calculated as z — 0.60 ± 0.11, making 
MMJ 04543 1.35 a possible cluster member. 

If MMJ 045421.55 and MMJ 045431.35 are cluster members 
at z = 0.54 SED-fitting suggests they each have Lfir ~ 5 x 
10 11 !/© and thus SFR~ 50 M yr~ 1 . We also find that they 
are colder than typical z ~ 2 SMGs, with dust temperatures of 
T d = 15±4K(for/3 = 1.5),orT d = 30±5K(for/3 = 1.1), com- 
pared to Td ~ 40 K and /3 = 1.5 for archetypal SMGs. Such prop- 
erties are not unpreced ented - the spectroscopic survey of SMGs by 
IChapman et alj J20051) contained examples of sources with equiva- 
lent millimetre-to-radio flux ratios at z ~ 0.5, suggestive of galax- 
ies containing cold dust. Similarly, in the SC UBA Local Universe 
Galaxy Survey (SLUGS) lDunne et al.1 feOOd) surveyed local IRAS- 
bright galaxies with SCUBA at 850 [im and utilised the combined 
IRAS and SCUBA photometry for SED fitting. This sample of lo- 
cal galaxies has T d = 35.6 ± 4.9 K, and (5 = 1.3 ± 0.2. In- 
deed, cold low-redshift galaxies are easier to detect at 850 \im or 
1.1-mm than their hotter counterparts. This is because colder dust 
produces emission which peaks at longer wavelengths than hot 
dust. Therefore, we do not find it unreasonable that one or both 
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of MMJ 045421.55 and MMJ 045431.35 are cluster members at 
z = 0.54 with T d ~ 30 K and (3 = 1.1. 

If both MMJ 045421.55 and MMJ 045431.35 are members of 
MS0451-03, their combined SFR is ~ 100 M yr -1 - a sig- 
nificant fraction of the SF R of all the clus ter galaxies within 2 
Mpc (200± lOOM vr~ 1 : lGeach et alj2006l) MMJ 045431.35 lies 
~ 2.6 Mpc from the cluster centre (about half of the turnaround 
radius), but MMJ 045421.55 is much closer to the centre: ~ 1 Mpc 
in projection. Notably, both of these systems are most likely merg- 
ers. Although not conclusive, if they are both cluster members, this 
suggests that existing starbursts are not instantaneously suppressed 
as they are accreted into the cluster environment, and can survive 
until at least 1 Mpc from the cluster centre. Alternatively the galaxy 
pairs could have been accreted into the cluster - suggesting that ac- 
cret ed galaxies can reta in their gas reserves during infall into clus- 
ters dGeach et al]|2009t) . 

3.4 Cluster and field SEDs 

We can also investigate obscured star-formation of the general 
galaxy population in MS 0451—03 below the flux limit of our 
AzTEC map by stacking fluxes at the positions of known cluster 
galaxies. Opt ical galaxies i n MS 0451—03 were morphologically 
classified by iMoran et al.l d2007bl) using the scheme defined by 
I Abraham et alj d 19961) . which we group into early-types and late- 
types for this analysis. We also define a mid-IR sample of galaxies 
which are bright at 24 urn, based on the catalogue of iGeach et al.l 
(2006) with 5*24^™, ^ 200 pJy. To reduce contamination we con- 
sider only those galaxies with spectroscopic redshifts and to search 
for environmental dependencies we also examine spectroscopically 
identified field galaxies. The cluster members are required to have 
0.52 < z < 0.56 and the field population is outside this win- 
dow. The field samples have median redshifts of 0.58, 0.46, and 
0.28 with interquartile ranges of z — 0.33-0.82, 0.26-0.62, and 
0.20-0.51 for the mid-IR galaxies, late-types and early-types re- 
spectively. For each of these six samples we stack the MIPS 24 (j.m, 
AzTEC 1.1 mm and VLA 1.4GHz fluxes. We note that the both 
galaxy samples are optical magnitude limited due to the require- 
ment for a spectroscopic redshift. Therefore, the stacked SEDs may 
not be representative of the entire population, in particular the most 
obscured galaxies are likely to fall below the optical magnitude 
limit. However, we expect such selection effects to equally affect 
the cluster and field samples allowing us to compare populations 
between the two density regimes. 

Any galaxies within 9 arcsec (the radius of the AzTEC beam) 
of an AzTEC map pixel with S/N^ ±3.5 are removed prior to 
stacking the AzTEC map. The contribution from each AzTEC-faint 
galaxy is weighted by the inverse of the squared noise at that pixel 
to calculate the weighted mean 1.1 mm flux of each population. 
Similarly we calculate the clipped weighted mean radio flux of each 
sample using our 1.4 GHz VLA map, correcting for bandwidth 
smearing, and stack 24 \im image in the same way. The stacking 
parameters are given in Tableland the resulting SEDs are pre- 
sented in Fig. [8] 

We fit the mid-IR and la te-type galaxies with template 
SEDs from iDale & HeloJ d2002l) and ca lculate the corre spond- 
ing SFRs based on the far-IR luminosity dKennicuttl Il998) . Mid- 
IR cluster galaxies without bright AzTEC counterparts have 
SFR = 16.5 ± 0.6 Mg,yr -1 which is consistent the SFR estimate 
of 17.0 ± 0.7 M yr _ for the mid-IR field population. In contrast, 
the late-type galaxies in the field have SFR = 3.0 ± 0.3 M yr _1 , 
compared to SFR = 2.3 ± 0.1 M yr _1 in the cluster. Both clus- 



ter and field early- type populations are undetected at 1.1 mm and 
24 p_m and SED fitting to the 24 p.m limits suggest, on average, 
SFR< 0.06 Mgyr -1 . 

Since our samples are flux limited, the different redshift dis- 
tributions will affect the derived SFR. The extent of this effect on 
our results is tested by calculating an observed SFR for each sam- 
ple by fitting templates with S2&\im = 200 [jJy at the redshift of 
the galaxies included in the stacks. Due to the generally higher red- 
shift of the field sample, we find that identical mid-IR populations 
would be observed with SFR 0.74 ± 0.16 times lower in the clus- 
ter than the field. In fact the stacked cluster population is observed 
with SFR 0.97 ± 0.05 times lower than the field population, there- 
fore, to within ~ la the mid-IR samples are consistent. However, 
the equivalent test for late-type galaxies, where the field sample is, 
on average, lower redshift, suggests that, due to the different lumi- 
nosity distance, identical populations in our analysis would appear 
2.6 times more active in the cluster than in the field. In fact, the 
cluster late-type galaxies are less active than the field. Therefore, 
the cluster environment is suppressing star-formation activity in the 
late- type population. 

IGeach et all d2006l) estimated the total SFR within 2 Mpc of 
the core of MS 0451-03 as 200 ± 100 M©yr _1 , based on con- 
verting 24 p.m fluxes of colour- selected 24-p.m det ected galaxies 
to tota l IR luminosities with the average SED from Dale & Heloj 
(2002). Our analysis presented in this paper includes radio and 1.1 
mm data enabling us to better characterise the cluster populations, 
and by stacking we can probe a fainter population. Within 2 Mpc 
of cluster core we calculate SFR> 315 ± 50 MQyr -1 from the 
spectroscopically confirmed mid-IR and late-type galaxies. To this 
we can then add 50 MQyr -1 for MMJ 045421.55, the potential 
ULIRG within 2 Mpc of the cluster core. 



4 CONCLUSIONS 

In this study we have investigated the dust-obscured star-forming 
population of the galaxy cluster MS 0451—03; we utilise 1.1 -mm 
observations to study obscured star-formation in MS 0451—03. We 
present a a ~ 1.1 mjy AzTEC map of the central 0.1 deg 2 of 
MS 045 1-03, within which 36 sources are detected at S/N> 3.5. 
We use radio, 24 (j.m, IRAC and SMA observations to precisely 
locate 18 of these SMGs. 

We calculate the reliability of photometric redshifts for SMGs 
and find that they are able to remove the bulk of background con- 
tamination, allowing us to isolate potential cluster members us- 
ing our optical, near- and mid-IR photometry. Based on these red- 
shifts we find two SMGs which are possible cluster members: 
MMJ 045421.55 and MMJ 045431.35. These systems are both re- 
solved into close pairs of galaxies by our ground-based and HST 
imaging, suggesting that interactions have triggered their starbursts. 
If they are cluster members both of these SMGs contain cold dust 
with T d ~ 30 K , for (3 = 1.1 (similar to SLUGS galaxies; 
iDunne et alfeood) and h ave SFRs of ~ 50 M Q yr -1 each. 

IGeach et aiT d2006 ) compared obscured activity, based on 
24|im emission, in MS 0451-03 and CI 0024+16 at z = 0.39, and 
found that CI 0024+16 has SFR ~ 5 times that of MS 0451-03 
within 2 Mpc of the clusters centres. They show that, taking into 
account the slightly different redshifts and masses of these two 
structures, MS 0451—03 is underactive and CI 0024+16 overactive 
at 24p.m. Therefore, it is likely the other galaxy clusters, including 
CI 0024+ 16, could contain ULIRGs in significantly larger numbers 
than we find in MS 0451-03. 
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Figure 8. Clipped weighted-average SEDs of the 2A\xm, late- and early-type populations of MS 0451—03 and the field. In the case of non-detections the 3<r 
noise limit is marked by arrows. The field populations are calibrated to their mean redshifts of 0.58 , 0.46, and 0.39 for the mid -IR galaxies, late- and early- 
types and respectively. We also show the best-fitting SEDs and corresponding SFRs I Kennicutt 1998) from Dale & Helou 1 2002) for the mid-IR and la t e-type 
galaxies; since the early-type populations are only detected at 1.4 GHz we cannot select one best-fitting SED and instead display all the Dale & Helot] (20021) 
SEDs on this panel. Although the observed populations have similar activity levels, in the late-type galaxies this is due to the flux limits of our observations 
and the different luminosity distances between the samples. If the cluster and field late-type populations were equivalent, the cluster should have SFR ~ 3 
times higher than observed. We deduce that the cluster environment has caused a reduction in the SFR of this population, and that these galaxies are probably 
in the process of transformation onto the red sequence. 



Table 5. Parameters and results of stacking on known cluster and field galaxies in the 24 \m\, 1.1 mm and 1.4 GHz images. SFRs are based on 
SED fits, except for the early-type galaxies, which are based purely on 24 \xm detection limits. Non-detections are represented by 3<r limits. 
Since the samples are flux limited and have various median redshifts we also present the expected observed SFR of the field mid-IR and late- 
type samples were they to have z mel a an = 0.54 ( i]3.4t . Differences between this value and the observed cluster SFR for the late-type galaxies 
suggests that the cluster late-type population is less active than the field late-type population. 



Population 


N 24u .ni 


o 6 

(MJy) 


Nl.lmm" 


a 6 

Jl.linm 

(MJy) 


Nl.4GHz a 


Q b 
01,4GHz 

(MJy) 


* median 


SFR 

(M yr -1 ) 


SFR (z = 0.54) 
(M Q yr -1 ) 


Cluster mid-IR galaxies 
Field mid-IR galaxies 


14 

84 


222 ± 18 
331 ±21 


14 

94 


< 970 
340 ± 130 


15 

89 


41.1 ± 1.6 
26.8 ±1.1 


0.54 
0.58 


16.5 ±0.6 
17.0 ±0.7 


16.5 ±0.6 

12.6 ±2.8 


Cluster late-type 
Field late-type 


167 
315 


38.4 ±7.7 
34.4 ±7.1 


131 

309 


< 340 

< 230 


162 
355 


5.2 ±0.6 
8.0 ±0.6 


0.54 
0.46 


2.3 ±0.1 
3.0 ±0.3 


2.3 ±0.1 
7.8 ± 1.1 


Cluster early-type 
Field early-type 


148 
35 


< 1.0 

< 2.1 


106 

32 


< 390 

< 530 


144 
40 


3.2 ±0.4 
10.9 ±0.9 


0.54 
0.28 


< 0.07 

< 0.03 


< 0.07 

< 0.19 



a The number of galaxies used in the stacking - excluding those close to the edge, off-image, or near detected sources 
^ Clipped weighted-average flux 



To further investigate the obscured star-forming population 
which lies below the limit of our AzTEC observations we create 
composite SEDs of spectroscopically confirmed mid-IR, early- and 
late-type cluster members and compare them to the corresponding 
field populations. As expected we find that both early-type popula- 
tions are undetected at both 24 \im and 1 . 1 mm and so are unlikely 
to be actively forming large numbers of stars (< 0.1 MQyr -1 ). 
The 24 Mm galaxies are significantly more active than the morpho- 
logically classified late-types with SFRs ~ 15 MQyr -1 versus 
~ 3 MQyr -1 on average. We find that the star-formation activity 
in the cluster late-type population, compared to a redshift-matched 
field population is quenched and ~ 3 times lower than expected. 
Mid-IR galaxies do not show this trend suggesting the more in- 
tense activity in these systems is more robust to environmental in- 
fluences. We find that the total SFR > 315 ± 50 M Q yr -1 in the 
central 2 Mpc of MS 0451-03. However, if MMJ 045421.55 is a 



cluster member it has SFR ~ 50 Mgyr and lies 1 Mpc from 
the cluster centre, taking the total SFR within 2 Mpc to > 360 

M yr -1 . 
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APPENDIX A: NOTES ON INDIVIDUAL SOURCES 

In Fig. lAll we present images centered on each AzTEC source, in 
order of decreasing S/N. We show Subaru and IRAC colour, and 
HST images, with radio and 24 |j.m contours. Identified galaxies 
are labelled and discussed below. 

1. MMJ 045438.96: SMA observations of this galaxy confirm the 
identified radio and mid-IR counterparts. The corresponding red 
optical and IRAC source at %,hot = 3.42^Q gg, has a disturbed 
morphology in the HST image. 

2. MMJ 045447.55: This galaxy has been located with the SMA 
but does not have any radio or 24 |i.m counterparts, or optical or 
IRAC galaxies at the SMA position. Therefore, MMJ 0455447.57 
could be one of a population of very distant SMGs, or be much 
cooler and more obscured than typical SMGs. 

3. MMJ 045433.57: SMA observations confirm the identified 
radio and mid-IR counterparts, corresponding to an optically faint 
but near-IR bright galaxy with z p hot = 2.55^g'gg . The HST image 
shows three components to this galaxy such that a merger is the 
likely cause of the ULIRG. 

4. MMJ 045431.56: We identify a robust extended radio counter- 
part 9.5 arcsec from the 1.1 mm centroid, which also corresponds 
to faint mid-IR emission; we consider the corresponding red galaxy 
at 2 p hot = 0.86^q'q3 to be the likely optical counterpart. The HST 
image of the region shows a galaxy with a bright core and a low 
surface-brightness tail or edge on spiral arm. However, the picture 
of this SMG is complicated by the presence of a statistically 
tentative mid-IR counterpart 9.1 arcsec from the 1.1 mm centroid 
which is unassociated with the aforementioned radio emission. 
The 24 [im position is coincident with an extended red IRAC 
galaxy but no optical sources. Further complexity arises from the 
SMA observation; MMJ 04543 1.56 was targeted but not formally 
detected. However, the SMA data does show two potential sources, 



close to an extended ~ 3<r peak in the radio map, ~ 3.7 arcsec 
from the AzTEC centroid and unassociated with both the radio 
and 24 \im counterparts previously discussed. There is also a red 
IRAC source at this position, and the HST image shows two optical 
galaxies. However, in the IRAC and ground-based optical data the 
counterpart is blended with a nearby elliptical galaxy. 

5. MMJ 045421.55: We find reliable coincident radio and 24 ^.m 
counterparts, corresponding to two optical galaxies separated 
by 2.6 arcsec, which are possible interacting cluster members 
(Zphot = 0.50to'?o)- The system is discussed in detail in ^3.3.31 

6. MMJ 045417.49: Both mid-IR and strong radio counterparts 
are detected although the presence of a nearby, bright star prevents 
multi-wavelength optical and near-IR study. Although the stellar 
halo is less extended in the HST image no counterpart is visible. 

7. MMJ 045413.35: The radio and 24p.m identification agrees 
with a 2-3cr peak in the 890 \im 'dirty' SMA map. Unfortunately 
this source lies outside of the IRAC and much of the ground-based 
optical coverage, although a galaxy is detected in the {/-band and 
HST imaging. 

8. MMJ 045412 .72: As potentiall y part of the gravitationally 
lensed arc from lBorvs etaU J2004bh this source is discussed fully 
in Wardlow et al. (2009b). O ur optical imaging i s too shallow to 
detect the ERO discussed in iBorvs et alj J2004bh . but radio and 
extended IRAC emission betrays the likely counterpart. 

9. MMJ 045345.31: There are two secure radio and one 24 [im 
counterparts. One of the radio galaxies corresponds to the 
24 [i.m position. The IRAC source has a red component at the 
mid-IR/radio position, but it blended with the emission from an 
optically saturated star making further conclusions difficult. It is 
unclear whether the 24 \im and radio flux is from the star or a 
background SMG counterpart which may be the cause of the red 
component of the IRAC emission. 

10. MMJ 045407.14: We find coincident reliable radio and 24 ^.m 
counterparts which are located in a red region at the edge of a 
bright, resolved galaxy which appears disturbed in the HST image 
and has z p hot — 0.351q'q3. The morphology of the galaxy sug- 
gests that a merger has triggered a dusty starburst region, causing 
the millimetre emission. The large size and optical brightness 
support the low redshift of this galaxy, although if the photometric 
redshift errors are underestimated it could potentially be a cluster 
member. Similarly, it is also possible that a background galaxy 
is the source of the millimetre emission. The detection of CO 
emission lines is required to confirm either scenario. 

12. MMJ 045426.76: There are no radio or mid-IR counterparts 
identified. Howe ver, an IRAC ga laxy has colours suggestive of 
SMG emission dYun et al] 1200 8h . yielding a counterpart with 
z phot = 0.89 ± 0.10. 

15. MMJ 045328.86: A 24 iim counterpart is detected 3.8 arcsec 
from the AzTEC position, corresponding to an IRAC source, 
but no detectable optical flux. The IRAC photometry yields 
Zphot = l-87lQ g5. There is coincident faint radio emission 
(~ 3cr) at the position of the 24 |j.m counterpart, increasing the 
likelihood that this is the correct identification. 

17. MMJ 045431.35: A radio counterpart 4.8 arcsec from the 
AzTEC position coincides with an pair of merging galaxies which 
are resolved in the HST image and have z p h t = 0.60 ± 0.11; this 
possible cluster ULIRG is discussed further in ^3.3.31 

18. MMJ 045411.57: 24-p.m emission betrays the source of the 
millimetre emission as a galaxy with z p hot = 0.74^o'?i, which 
is resolved in the HST image as an interacting pair, separated by 
< 1 arcsec. There is faint radio emission (~ 3<r) ~ 2.3 arcsec 
to the northeast of the 24 ixm identification, coincident with a red 
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Figure Al. 25 X 25 arcsec images centred on each AzTEC galaxy are shown in the left-hand and middle panels; north is up and east is to the left. In 
the left-hand panel we show Subaru data: BVR colour images with solid radio contours at -3 (dotted), 3, 4, 5, 6, 7, 8, 9, lOxcr are overlayed. The middle 
panel contains true-colour images containing the IRAC 3.6 + 4.5 (blue), 5.8 (green) and 8.0 |jm (red) data; dashed contours present 24 \xn\ flux at 5, 10, 
20, 30, 40, 50, 100, 200, ...1000X Ljjy per pixel and the right-hand panel contains 20 X 20 arcsec //STF814W cutouts centered at the AzTEC position. In 
all images circles are centred on the AzTEC positions and have 20 arcsec diameter, corresponding to our radio and mid-IR search radius. 24|am and radio 
counterparts are highlighted with diamonds and squares respectively; dotted symbols represent tentative counterparts (with 0.05 $J P SC 0.10) and the sizes 
of the symbols are r epresentative of the typical astrometric errors of these data. We mark the positions of SMA detected sources with stars and X-ray sources 
lMoln aretalj2002l) with crosses. 
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Figure Al - continued 



IRAC source. We conclude that the interaction between the optical 
galaxy pair is producing all of the emission. However, it is possible 
that the radio and IRAC source is unrelated to the optical emission, 
in this case we can only constrain the redshift to z p hot > 1.14. 
26. MM J 045349.69: A faint optical galaxy with z ph ot = 
1.33+ 

aio corresponds to the robust mid-IR counterpart. The 



galaxy is undetected by HST and the region is not covered by our 
IRAC mosaic. 

27. MMJ 045421.17: The mid-IR counterpart is coincident 
with both faint radio emission and a faint red galaxy at 
Zphot = 1-801q28, which the HST imaging resolves into a 
bright nuclear region and some extended emission. 
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Figure Al - continued 



28. MMJ 045345.06: The galaxy is tentatively identified through 
its mid-IR emission, which is coincident with faint radio flux; 
our IRAC mosaic does not cover the counterpart, but optical 
photometry yields z p hot = l-SQ^o'sl- 

29. MMJ 045442.54: Although there are no radio or mid-IR 
counterparts MMJ 045442.54 is identified from its IRAC colours 
and the photometry of the corresponding faint galaxy yields 
zphot = 1.09 ±0.15. 
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